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Abstract 

We compute the probabihty distribution for bubble cohisions in an inflating false vac- 
uum which decays by bubble nucleation. Our analysis generalizes previous work of 
Guth, Garriga, and Vilenkin to the case of general cosmological evolution inside the 
bubble, and takes into account the dynamics of the domain walls that form between 
the colliding bubbles. We find that incorporating these effects changes the results dra- 
matically: the total expected number of bubble collisions in the past lightcone of a 
typical observer is ~ 7 Vf /Vi , where 7 is the fastest decay rate of the false vacuum, 
Vf is its vacuum energy, and Vi is the vacuum energy during inflation inside the bubble. 
This number can be large in realistic models without tuning. In addition, we calculate 
the angular position and size distribution of the collisions on the cosmic microwave 
background sky, and demonstrate that the number of bubbles of observable angular 
size is Nls ~ V^^-^j where fifc is the curvature contribution to the total density at 
the time of observation. The distribution is almost exactly isotropic. 
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Figure 1: The Persistence of Memory, Salvador DaK, 1931 




Figure 2: The Disintegration of the Persistence of Memory, Salvador DaK, c. 1952-1954. 
"My father today is Dr. Heisenberg" -Salvador Dali, Anti-Matter Manifesto, 1958 
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1 Introduction 



In a model like the string theory landscape [1, 2] one may expect our observable universe 
to be contained inside a bubble which is surrounded by an eternally inflating false vacuum. 
Quantum tunneling and thermal effects allow the false vacuum to decay, producing regions 
of lower energy — like the one we inhabit — which expand and may themselves decay in turn. 
If the decay rate of the false vacuum in Hubble units is less than a particular order one 
number, this process never terminates. The false vacuum produces new volume through 
exponential expansion faster than it loses volume due to decays. 

Nevertheless any bubble will collide with many others, forming a cluster that grows 
without bound as time passes, as first described in [3]. We would like to ask whether 
observers like ourselves can expect to have a bubble collision in our past. If the answer is 
yes, then we would like to know the expected number of collisions, their distribution in size 
and direction, and their observable consequences. 

This distribution of vacuum bubble collisions was analyzed recently in a very interesting 
paper [4] by Guth, Garriga, and Vilenkin (GGV) titled "Eternal inflation, bubble collisions, 
and the persistence of memory." The analysis of GGV was done in the approximation that 
the cosmological constant inside the observation bubble is the same as the cosmological 
constant of the false vacuum. Also, GGV assumed that observers cannot form anywhere in 
the future lightcone of coUisions. 

With these assumptions, GGV computed the probability distribution for bubble collisions 
which may become visible in the future. They found that a typical observer inside such 
a bubble sees an anisotropic distribution of collisions. The anisotropy is left over from 
the initial conditions, and does not go away even in the limit that the observation bubble 
nucleates an inflnitely long time after the initial conditions. They called this Angering effect 
of the initial conditions "the persistence of memory." 

However they also found that the overall rate at which an observer sees collisions (and 
the probability of observing any collisions) is proportional to the dimensionless decay rate 7 
of the false vacuum, a quantity that is typically exponentially small. Moreover they found 
that the typical distance to the nearest collision is large, of order — In 7 in units of curvature 
radius. 

In this note we consider the case relevant for observation, where the interior of the 
observation bubble is similar to our universe. We assume that due to the dynamics of the 
domain walls, collisions disrupt structure formation in only part of their future lightcone (as 
described in more detail in section 2.1). We flnd that these changes have a major impact on 
the result. The ratio of the vacuum energy in the false vacuum to the inflationary vacuum 
energy, Vj/V^, which was set to one in the GGV approximation, can be large in realistic 
scenarios. Among other things, this ratio turns out to divide the anisotropic term in the 
distribution and multiply the expected number of bubbles in the past lightcone of a late- 
time observer. Incorporating the effects of domain walls renders the distribution smooth and 
finite even for observers at infinite "boost" . 

A heuristic argument for the distribution of bubbles in our past lightcone is as follows. 
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In the inflating false vacuum there is an event horizon with radius oi Rf — ^ (where 
Hf ~ \/~^f is the false vacuum Hubble constant) surrounding every point. This horizon 
shields the bubble from interacting with any event (such as another bubble's nucleation) that 
occurs more than a distance Hj^ away from its walls. Hence the bubble wall is surrounded at 
all times by a shell Hj^ thick inside of which any bubble nucleation leads to a collision. As the 
bubble grows the volume of this region expands like its surface area, and at the same time the 
relative size of the late-appearing bubbles shrinks. This leads to a size distribution of bubble 
collisions which, as a function of angular radius on the fiducial bubble, increases rapidly with 
decreasing size. Further, as bubble grows the initial conditions become unimportant, so the 
anisotropy is insignificant for small bubbles. 

To compute the total number of bubbles, we need to know the 4- volume in the false 
vacuum available to nucleate bubbles that will collide with our own. The surface area of the 
bubble wall inside the past lightcone of late-time observers like us is of order H^^, where 
Hi is the Hubble constant during slow-roll inflation inside the bubble. As mentioned above, 
the distance from the bubble wall to the event horizon is Hj^. To get the 4- volume, we just 
need to know the amount of time available for creating collision bubbles, which turns out to 
also be given by Hj^. 

As a result the 4-volume of the false vacuum available for bubble nucleations scales as 
Hj^H~^. If r is the decay rate per 4-volume, the total number of bubbles in our past 
lightcone is ~ TH^^HJ^ Defining the dimensionless decay rate 7 = THJ*^ we have 

iv~7(iy. (1.1) 

It is crucial in deriving this answer that observers can form in a reasonable fraction of the 
future lightcone of the collision. Whether this is true depends primarily on the dynamics of 
the domain wall that is formed in the collision. The domain wall will have constant proper 
acceleration and approach the speed of light. Depending on parameters it may accelerate 
into our bubble or away from our bubble. 

If the domain wall accelerates into our bubble, then most of the future lightcone of a 
collision is uninhabitable because it is eaten by the collision bubble, and the probability 
of observing a collision will be small, as in GGV. Further, such an observation would be 
followed nearly instantaneously by a fatal collision with the domain wall. 

We choose to analyze the more optimistic case where the domain wall formed in the 
collision accelerates away from our bubble. As we describe in section 2.1, it is reasonable 
to expect that in this case an order one fraction of the future hghtcone of the coUision is 
inhabitable, and the estimate (1.1) holds. 

Finally, one could also analyze the case where our bubble collides with identical bubbles. 
In this case no domain wall forms after the collision, and again we expect an order one 
fraction of its future lightcone to be inhabitable. We will focus on the case where the collision 
produces a domain wall, but we expect that our analysis carries over almost unchanged to 
the collision of our bubble with identical bubbles. 
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1.1 Results 



Our primary result is that the expected number of bubble collisions in our past is 



Here 7 ~ e^'^ is the dimensionless decay rate of the false vacuum to its fastest decay product. 
The inflationary Hubble scale Hi is very poorly constrained, but assuming the vacuum energy 
in the false vacuum is of order the Planck scale, observational bounds on tensor modes imply 
that {Hf/HiY > 10^^, and it can be larger than 10^° in models of low scale inflation. Also, 
the decay rate is bounded below by the recurrence time of the false vacuum de Sitter space, 
7 > exp(—Sf) with Sf ~ {Mp/ Hf)^ . If the vacuum energy of the false vacuum is near 
the Planck scale the decay rate cannot be too small. So even though bubble nucleation 
is exponentially suppressed, the total number of bubbles in the past lightcone of a typical 
observer can be large with no fine-tuning. 

We derive the distribution of bubble collisions in the backward lightcone of a single 
observer as a function of the angular coordinates 6 and of the center of the collision 
bubble, a parameter X that controls the center of mass energy of the collision, and the 
conformal time rjy when the collision first becomes visible (Eq. 2.9). 

We find that the "persistence of memory" effect discovered by GGV is significantly modi- 
fied when the vacuum energy during inflation is allowed to be different from the cosmological 
constant of the false vacuum. The memory of the initial conditions, which manifests itself as 
an anisotropy in the angular distribution, disappears in the limit that the bubble has small 
cosmological constant inside (Eqs. 2.9 and C.23) and the late-time distribution becomes 
precisely boost and rotationally invariant around the bubble's nucleation point (Eq. B.39). 

Finally, we derive the distribution of collision lightcones on the last scattering surface as a 
function of their position and angular size on the CMB sky (Eq. 4.25). Bubbles that collide 
with ours very early have lightcones that completely cover the part of the last scattering sur- 
face we can see, and these are probably not observable in the limit that the size of the region 
they cover is much greater than the observable part. But bubble collisions that occurred 
somewhat later have lightcones that bisect the observable part of the last scattering surface 
(and hence the CMB sky), and these can directly affect the cosmic microwave background 
temperature map and are of observational interest. The total number of such bubbles is 



7Vi5 ~ 4CoiVV^^fe = — Co7 J V^k. (1.3) 

where ^ 10/3 is a constant the depends upon the cosmological evolution inside the bubble 
after reheating (see Appendix D for details) and flk is the curvature contribution to fltot- 

1.2 Related work 

Our work is closely related to the analysis of Aguirre, Johnson, and collaborators [5, 6, 7], 
with the crucial difference that we treat the effects of bubble collisions on observer formation 
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in a different (and in our opinion more realistic) way. Our results are in rough agreement 

with [5] for the number and distribution of "small" bubbles, where "small" here means bub- 
bles that are smaller than the observation bubble in the observer's reference frame. The 
number of such bubbles was computed in [5] and their isotropy was noted by [5] and [8]. 
However because of our assumptions about the disruptive effects of collisions on observer for- 
mation, our results for the distribution of large bubbles and our conclusions for the potential 
observability of collisions differ sharply from [5, 6, 7]. 

Ref. [5] focused on the case where bubble collisions have approximately no effect on 
observer formation, and found that typical observers see an arbitrarily large number of very 
large, anisotropically distributed bubbles. The infinity is removed here by our treatment of 
the domain walls. Ref. [7] did a numerical analysis of a single collision and found that the 
region far to the future of the resulting domain wall has an approximate SO{3, 1) invariance 
restored (similar to the case of a single bubble without a collision), which may remain if no 
other decays or collisions take place to perturb it. Based on this and on the conclusions of 
GGV, [7] concluded that the fraction of observers that can see a collision is very small. How- 
ever, even under this set of assumptions computing this fraction unambiguously is difficult, 
as it requires comparing two infinite volume regions. 

After this work was completed and while the manuscript was in preparation we received 
Ref. [9], which concluded that the fraction of observers that see collisions can be C(l), but 
performed the analysis in the approximation that the domain walls between bubbles do not 
accelerate and that the observation and collision bubbles are identical and dominated by 
vacuum energy inside. In this work we will take the acceleration of the domain walls into 
account and allow for a fully realistic cosmology inside the observation bubble. 

The motion of the post-collision domain walls was studied in [10] and [11] in special cases, 
and then more generally in [12] and [6] . 

Determining the observational consequences for observers in the space to the future of a 
bubble collision is a rich problem that has been addressed by [5, 6, 7] and by [12, 13]. In 
particular [13] studied in detail the effects on the cosmic microwave background of collision 
lightcones that bisect the observable part of the last scattering surface. 

2 Probability distribution for bubble collisions 

We want to compute the probability distribution for bubbles that collide with our own. 
Before including other bubbles, the exterior of our bubble is in the false vacuum. We work 
in the approximation that the domain wall of our bubble is lightlike, and that the metric 
outside the lightcone is undisturbed de Sitter space with cosmological constant Hj. This 
approximation does not have a large effect on our results, and going away from it would 
introduce new model-dependent parameters. 

A convenient coordinate system that covers the region where collision bubbles can nucle- 
ate is 



ds^ = H 



f cosh^A" 
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where dS^ is the metric of a unit 2 + 1 dimensional de Sitter space and Hf is the Hubble 

constant of the false vacuum. Boosts that leave the observation bubble fixed act as the full 
symmetry group of the 2 + 1 de Sitter space, while leaving X fixed. Therefore collisions with 
different X are physically different, while any two collisions with the same X can be boosted 
into each other by acting with the symmetries. 

Choosing coordinates on the de Sitter space, the metric is 

ds^ = Hf ^^{dX^ - dr^ + cosh^ rdfll) (2.2) 

cosh X 

where X and r run from —oo to oo. These coordinates cover the region outside our bubble 
that is capable of nucleating bubbles that will collide with our own. The nucleation point 
of our bubble is at A" = — oo and r finite, while the domain wall of our bubble is the null 
surface X = —oo with r > 0. The null surface X — oo bounds the region where nucleated 
bubbles will collide with our bubble. 

We are interested in observing bubble collisions from inside the bubble, while the co- 
ordinate system above covers the outside of the bubble. For the moment we leave the 
cosmological history general inside the bubble. Before considering the effects of collisions, 
the interior geometry has SO (3, 1) symmetry. The metric can be written 

ds'^ = a^{'q){-drf + dp^ + sinh^ pdVtl), (2.3) 

where without loss of generality we focus on an observer at p = 0. This metric describes 
a homogeneous and isotropic Robertson- Walker cosmology. An unperturbed observation 
bubble will inflate and reheat along slices of constant r/, and comoving observers move on 
trajectories with fixed spatial coordinates^. 

This metric is smooth all the way down to the "big bang" r] = —oo, even though a{r]) ^ 
there. The reason for this non-singular "big bang" is that the instanton boundary conditions 
require [14] 

a{r]) ^ -^e" + Oie^") as 77 ^ -oo (2.4) 

where the multiplicative factor 2/Hf is arbitrary and has been chosen for convenience. For 
example if one follows the past lightcone p -|- 77 = 770 of an observer at time 770, the radius of 
the bubble hghtcone is 0(77) sinh(p) — > e'^'^^/Hf — e^°/Hf, which is finite as 77 — > —00 and 
p — > 00. 

Consider a bubble that nucleates at a point {X ,t,9,4>). The coUision occurs along a 
spacelike hyperboloid that preserves S0{2, 1) symmetry, as shown in figure 3. The bubble 
will appear as some type of disturbance on the observer's sky. The angular location of the 
center of the disturbance is given by the angular location {6, 0) of the bubble nucleation. 
We would like to compute the time at which the collision first enters the backward lightcone 
of the observer at p = 0. Ingoing radial null rays satisfy X + r = c outside the bubble and 

^Notc that we are using conventions where the scale factor a has units of length, and the spatial curvature 
is fc = —1. 
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Figure 3: Two views of tlie collision of two bubbles in spacetime. The red cone is the observation 
bubble lightcone, the green cone is the collision bubble lightcone, the straight black line is the 
worldline of the observer, and the hyperbolic black line is the surface along which the bubbles 
collide. The angular size of the collision as seen by the observer at some time is determined by the 
intersection of the observer's past lightcone (in blue) with the observation bubble wall, or with the 
surface of last scattering. 



(with our choice of coordinates) t] + p = c inside the bubble, where c is an arbitrary constant. 
So the future lightcone of a bubble nucleated at {X, r) first reaches p = at the conformal 
time 

r], = X + T. (2.5) 

A given collision is described by four parameters: two angular coordinates {6, 0) give 
the angular location of the center of the collision, rjy is the time when the collision first 
becomes visible, and X sets the intrinsic properties of the collision. For example, the radius 
of curvature of the collision hyperboloid is Vc = exp(A'). 

The probability to nucleate a bubble in an infinitesimal region is proportional to the 
4-volume of that region, 

dN = -fHjdVi, (2.6) 

where 7 is the dimensionless nucleation rate. The infinitesimal 4-volume in our coordinates 
(2.2) is 

dVi = H7^ ^ cosh^ TdTd'^^2- (2.7) 

^ cosh X 

We can exchange the coordinate r for the conformal time rj^ at which the collision enters 
the backward lightcone by using (2.5). So the naive probability distribution for collisions is 

dN = ^ cosh^(^^--y) ^^^^^^2^^ (2.8) 
cosh X 
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Again naively, the total number of collisions visible to an observer at conformal time rjo is 

given by integrating the above equation over the entire range of parameters, subject to the 
constraint rj^ < rjQ. However this integral diverges at rj^ — oo, corresponding to an infinite 
probability to have nucleated a bubble in the arbitrarily far past. 

There are several reasons the above analysis is naive, which is fortunate since the answer 
is infinite. We discuss various corrections throughout the rest of the paper. 

2.1 Disruptive effects of collisions 

The most important correction is that bubble collisions will disrupt structure formation in 
some part of their future lightcone. Since we are interested in making a prediction for our 
observations, we need to compute the distribution of bubble collisions consistent with the 
presence of observers. We show in this subsection that this requirement regulates the infinity 
of the naive analysis and allows for a well-defined prediction for the distribution of bubbles.. 

The collision can disrupt structure formation throughout its future lightcone, but it is 
reasonable to expect that in some regions the bubble collision has a small effect, while in 
other regions no structure forms. We focus on the collision of our bubble with some other 
type of bubble, so that a domain wall forms in the collision. Vacuum energies and domain 
wall tensions are generically Planckian in the landscape, so we assume that both the false 
vacuum and the collision bubble have Planckian vacuum energy, and that the domain walls 
in the problem also have Planckian tension. The vacuum energy in the collision bubble could 
be negative or positive. 

We make the approximation that structure formation occurs unchanged in the observation 
bubble, except that the domain wall with the collision bubble removes regions that otherwise 
would have formed structure. Because the vacuum energy in the collision bubble is Planckian, 
no observers can form in the collision bubble. 

The key question is regarding the dynamics of the domain wall. Because we have assumed 
that the vacuum energy of the collision bubble and the tension of the domain wall are 
Planckian, the acceleration of the domain wall is generically Planckian. But the domain wall 
may accelerate into the observation bubble, or away from the observation bubble. In either 
case the domain wall undergoes constant proper acceleration and asymptotically appraches 
the speed of light. 

If the domain wall accelerates into the observation bubble, then nearly all of the future 
hghtcone of the collision is behind the domain wall, and since we have assumed that observers 
can only form in the observation bubble there are essentially zero observers in the future 
lightcone of the collision. 

If the domain wall accelerates away from the observation bubble, then some observers will 
form in the future lightcone of the collision. We focus on this case. If the vacuum energy of 
the collision bubble is positive, then it is guaranteed that the domain wall accelerates away 
from the observation bubble. If the vacuum energy of the collision bubble is negative, then 
the acceleration of the domain wall depends on the tension; for sufficiently large tension the 
domain wall again accelerates away from the observation bubble. 
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The detailed dynamics of the domain wall are somewhat involved and model-dependent 
[10, 11, 12, 6], but with our assumptions Planckian acceleration is generic. The initial 
velocity of the domain wall after the collision is also model dependent, but typically it will 
take some time for the messy collision to settle down into a coherent domain wall; with 
our assumptions this time scale is also Planckian. In order to have some control over our 
calculation we assume that the energy scales are all somewhat less than Planckian, and 
instead set by the outside Hubble constant Hf.'^ 

Therefore, we adopt the following caricature of the domain wall motion: we assume the 
domain wall moves into the observation bubble at the speed of light for a time Hj^. Then 
it turns around and moves outward at the speed of light. While this approximation could 
certainly be improved upon in specific models, it captures the crucial physics of the problem. 

As discussed in the paragraph below Eq. 2.3, the metric has a coordinate singularity 
a — >■ as — > —oo. With our choice of coordinates 1] = corresponds to a time of order 
Hj^. Therefore we will approximate the domain wall trajectory as null and ingoing for 
77 < 0, null and outgoing for 77 > 0. Bubbles whose domain walls cover p = at any time 
during their evolution are assumed to disrupt reheating. Therefore, given our assumptions, 
we prohibit bubbles with 77^ < 0. 



dN = 7 \ 4 U xdr]^(fQ2 (2.9) 



2.2 Summary 

Now we can compute the distribution of bubbles that are in the backward lightcone of an 
observer at the conformal time r]o. It is 

cosh^(r;^ — X) 
cosh^ X 

< Vv<Vo (2.10) 

where the restriction < r^t, comes from requiring the observer formation is not disrupted. 
We can integrate out X to get the distribution as a function of conformal time, 

27 

dN = -^(1 + 2cosh2rjy)dri^d^n2 for rjy > (2.11) 

Integrating, the total number of collisions visible before some conformal time 770 is 

N{vo) = -j^ (sinh 27^0 + 7^0) ■ (2.12) 

There is an important subtle point that we have glossed over here. The divergence present 
in the naive distribution means that the probability of observer formation at p = is zero. 
Therefore, one could ask whether the above distribution is meaningful. The answer is that 
any choice of initial conditions will regulate the infinity so that the probability of observer 
formation is finite. Further, the actual choice of initial conditions will have little effect on the 
distril)ution. so that the net effect is simply that we can take the above equations seriously. 



^In fact, the typical energy scale associated with the vacuum energy, the tension of the wall, and its 
acceleration, is the geometric mean between the Planck scale and Hf. In the following we will ignore this 
point, for it does not lead to any substantial modification of our results. 
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3 Effect of the initial condition surface: how persistent is mem- 
ory? 

Because the false vacuum is only metastablc, it is not consistent to imagine that it exists 
eternally. To have a well-defined Cauchy problem, one needs to choose an initial condition 
on a spacehke or null surface. An obvious (but ad hoc) choice is to put the fields in a false- 
vacuum minimum on some slice of constant de Sitter time and then consider the forward 
time evolution. One of the surprising results of the work of [4] was that observers inside 
bubbles that form to the future of this initial condition surface always retain some memory 
of it no matter how far in the future they live, in the sense that the angular distribution 
of incoming bubbles remains significantly anisotropic except for the special single observer 
in the case of an isotropic initial condition surface who sees an isotropic distribution. We 
will see that incorporating a realistic cosmology inside the bubble in the analysis alters this 
conclusion significantly. 

For an initial condition we follow [4] in taking the surface to be i = — oo in the standard 
de Sitter fiat slicing. This choice has the advantage that it is infinitely far back in the past, 
but still regulates the divergence in Eq. (2.7). On the other hand, [4] found that with this 
choice only a single special observer co-moving inside the observation bubble will see an 
isotropic distribution, because generic co-moving observers are boosted with respect to the 
flat slicing. 

For this discussion, since we are interested in comparing different observers related by 
boosts around the nucleation point of the bubble, it will be helpful to use the embedding 
coordinates for de Sitter space in which all of the symmetries are manifest. 3 + 1 dimensional 
de Sitter space can be embedded in 4 -|- 1 dimensional Minkowski space as the surface 

X^ + R'^-T^^l (3.1) 

where R"^ — Y'^ + Z"^ + W"^ and we have set 

Hf — I'vci this section. (3.2) 

The initial condition surface in embedding coordinates is X -|- T = 0. 
The relation to the parameters of the bubble coUision is 

^ sinh r 

cosh X 

X ^ - tanh X 

coshr 
Y = — — - cos 6 
cosh A. 

coshr 

Z — — -— — sm cos 
cosh X 

coshr . „ . , 

W = — — sm 6 sm 

cosh 
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where we have chosen the observation bubble to nucleate at X — 1,T — Y — Z — W — 0. 
Recall that the coordinate r is related to our preferred parameter 1].^ hj 7]^ = t + X . 

With this choice for the initial condition surface, allowed collision bubbles nucleate in 
the region X + T > 0. Translating to the parameters of the bubble collision, this condition 
is 

smh. x 

- tanh X H — > , (3.3) 

cosh X 

which is equivalent to 

X <T . (3.4) 

This is the appropriate condition for the special "unboosted" observer (the observer that is 
comoving in the bubble and at rest in the flat slicing). To consider boosted observers, we can 
simply boost the initial condition surface. The boost should preserve the nucleation point 
of the observation bubble; consider without loss of generality a boost in the YT plane. The 
boosted initial condition surface is 



X + Tcosh^ + ysinh^ > (3.5) 
where ^ parameterizes the boost. In our coordinates this condition becomes 

sinh X < sinh r cosh ^ + cosh r sinh ^ cos 9 . (3.6) 

The 4-volume available for bubble nucleations is restricted by the above inequality. 

Now the initial condition surface is ambiguous, because we do not know our boost 
However, we will show that the expected number of bubbles that satisfy the other constraints 
but are prohibited by the initial condition surface is generally much less than one! Therefore, 
the unknown parameter ^ has a small effect on the final distribution. 

To show this, we will compute the total number of collisions in the past lightcone of an 
observer in a cosmology like ours that are allowed by the other constraints but prohibited by 
the initial condition surface. The answer for arbitrary ^ is easy to calculate but is slightly 
complex, so wc will just present the answer at zero and infinite boost. 

At zero boost the initial condition surface imposes X < t, so the expected number of 
bubbles that are allowed by the other constraints but prohibited by the initial conditions 
surface is 

AN{i = 0) = jJdX JdrJ d'n^^ (3.7) 
integrated over the region 

0<A' + T<77o X >T (3.8) 

Evaluating this integral gives AA^(^ = 0) = (47r7/3)(l + ln4) + 0(76"'"'). 

At infinite boost, the initial condition surface requires cos ^ + tanh r > 0. Performing the 
integral, the number of bubbles prohibited by the initial condition surface only is AA^(^ = 

00) = (47r7/3)(r7o + 1) + 0{-fe-^^"). 



11 



We will show in the next section that the conformal time 770 for an observer like us is 
approximately 

770 ~ log ^ (3.9) 

where is the Hubble constant during slow roll inflation. Since one expects 7 to be 
exponentially small, generically AA^(^) <S 1 V ^. 

4 Observability of Collisions 

In this section we will use this result to determine the total number of collisions that are 
potentially detectable by an observer in a realistic cosmology. 

4.1 Computation of the conformal time 

To justify our claim that r/o ~ log we need to relate 77 to the proper time of a comoving 
observer in the observation bubble. The conformal time 77 is related to the observer's proper 
time t hy dr] — dt/a{t). The constant of integration is fixed by demanding that 

as 77 ^ -00 (4.1) 

as mentioned in the discussion below Eq. (2.3). 

To proceed we will need an explicit form for a{t). Approximating slow-roll inflation by 
de Sitter space with Hubble constant iJj, a{t) — Hf^ sinh(ifji) for t < t^. During this time, 
the scale factor as a function of conformal time is 

a\ri) = „ ^ (4.2) 

smh (7/ - log ;^) 

where the additive constant in 7] is fixed by the convention (4.1). 

Defining the number N of efoldings by a^H = H^^e'^ , the conformal time at reheating is 

VRH = log ^ - e-^ ^ log ^ (4.3) 

for ttRH > H~^. 

After reheating, the behavior of the scale factor depends upon which type of matter or 
energy is dominating the expansion of the Universe at a particular epoch. It is approximately 
a{t) — H^^e^ {t/tRHY ■, where p = 1/2 or 2/3 for radiation or matter domination respectively. 
Eventually a will asymptote to t (if there were no dark energy) or to an exponential if there is 
a cosmological constant in the bubble. We consider in detail the realistic case with arbitrary 
matter, radiation, curvature, and cosmological constant content in Appendix D. Here, we 
illustrate the calculation for the simpler case where the scale factor is a power law of the 
cosmic time. 
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Taking a{t) — H- ^e^ {t/tuHY and integrating, we get 

t t. 



,W-,H«=^^^-^j (4.4) 
For times much later than reheating, and many efoldings of inflation, this is 

We would like to relate this to observable quantities. The "curvature contribution" to the 
total energy density is defined by 

We have a — pa/ 1 so the conformal time can be written 

77(t)=log^ + -^xAM^) (4.7) 

Finally, most of the time (and most of the conformal time) since reheating is during matter 
domination, so setting p = 2/3 we have 

rj{t)^\og^ + 2,/ih{tj (4.8) 

valid for times well into matter domination and before dark energy or curvature domination. 
More generally we find that the expression for the conformal time today is 

r?o~log-^ + CoV^^fc (4.9) 

where Co ~ 3.3 for allowed values of the cosmological parameters (see Appendix D). 

The conformal time grows logarithmically with time for the first efold of infiation, rapidly 
asymptoting to log Hj j Hi. After inflation it remains almost constant until late times, when — 
if the curvature becomes significant — it will begin to grow again. If our universe has a small 
positive vacuum energy and is nearly flat then curvature domination never occur and the 
effects of collisions weaken from the onset of cosmological constant domination onward, 
making the current cosmological era an opportune time to observe them [15]. 

The total number of collisions in a universe similar to ours at the present time will be 



87r7 , . , „ , 47r7 

Is fa —— (smh 2r]o + rjo) ^ —— 



M e^^°V^ + 21og^ + 2Coy^W 



(4.10) 



To compute the total number in our backward lightcone now, we use flk{to) -C 1 and also 
expand for 3> 1 to get 

(§^) ' ^1 + 2CoV^ + ...) (4.11) 
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where the ... denotes terms which are much smaller than one. (See appendix C of [6] for 
a more detailed computation of this quantity.) The V^H^ term is at most an order one 
multiplicative correction, so the bottom line is 

There is no reason to expect this number to be small. While one expects 7 <C 1, observational 
constraints require [Hj/HiY > 10^^ if Hj ~ Mpi. In models of low-scale inflation Hi can be 
very small; in such models the ratio {Hf/Hi)'^ can be 10^'' or more. Moreover, the relevant 
7 is the fastest decay channel of the false vacuum around the observation bubble. In models 
like the string theory landscape, there are an enormous number of decay channels available. 
7 < 1 is required for our analysis to be valid, but it is not clear how much smaller than one 
we should expect it to be. 

4.2 Distribution at Last Scattering 

One of the most promising possibilities for observing these collisions comes from the cosmic 
microwave background. In order for a collision to make an observable mark on the CMB, the 
collision bubble lightcone should be of reasonable size when it intersects the part of the last 
scattering volume we can observe. If it is smaller than a fraction of a degree it will be lost in 
the noise, and if it is much larger than the observable part of the reheating volume the signal 
will be degenerate with the dipole from the peculiar velocity of the earth in the CMB rest 
frame [13]. However if it has an angular size of order one, it can create potentially observable 
disks on the CMB inside of which the average temperature is affected in an angle-dependent 
way [13]. Therefore the collisions of most interest observationally are those with lightcones 
that bisect the observable part of the last scattering volume. 

We describe the bubbles first in terms of their comoving distance from p = at last 
scattering. Then we will convert this to their angular size. For a collision which is not yet 
in the backward lightcone of p = at last scattering, the shortest comoving distance to the 
lightcone is given by the part of the lightcone which travels radially, 

P + riLS = riv (4.13) 

so that the comoving distance is 

^iVLs) = Vv - Vls (4.14) 

For bubbles that are already in the backward lightcone of p = at last scattering, i is 
negative but its magnitude still gives the shortest comoving distance from p = to the 
collision lightcone. 

To summarize, we will trade in 77^, the time the future lightcone of the collision crosses 
p = 0, for i, the comoving distance from p = to the collision lightcone at last scattering. 
We allow i to be positive or negative, with positive £ corresponding to collisions that have 
not yet crossed p = at last scattering, and negative i corresponding to collision that have 
already crossed p = 0. 
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Let's rewrite our distribution in terms of i. We change variables in the distribution by 
substituting 

nv = ^^y]Ls (4.15) 

so the distribution is 

dN ^ ^ '''^'^^^f'-^^ dxm^n, (4.16) 

cosh X 

< i + VLS (4.17) 

where the lower equation arises from demanding that the collision does not disrupt structure 
formation at p = 0. 

We can integrate out X as before to get 

dN ^^[1 + 2 cosh{2{i + r)Ls))]did'^n2 (4.18) 

This gives the distribution of comoving distances to coUision lightcones at last scattering. If 
we were interested in the disribution at some other time rji, we could just replace tjls Vi 
in this formula. 



4.3 Effects on the Last Scattering Surface 

A key measure of the observability of a given collision is the the ratio pis/^, where p^s is 
the comoving distance to the last scattering surface and i is the comoving distance to the 
collision lightcone. Recall that i is negative for collision lightcones that have crossed the 
origin at last scattering and positive for coUision hghtcones that have not yet crossed the 
origin. 

Bubbles with i > pis are spacelike separated from the entire last scattering surface and 
have no effect. Bubbles with —pis < ^ < Pls affect only part of the last scattering surface 
and have the most striking signals. Bubbles with £ < —pLs influence the entire last scattering 
surface. They may lead to interesting effects, especially at long wavelengths. However, it 
is clear that as the distance to the lightcone gets much larger than the distance to the last 
scattering surface, I/ Pls ~^ —oo, the anisotropy induced by the collision must vanish. 

The comoving distance to the last scattering surface is related to the curvature, 

PLS ~ CoVMto), (4.19) 

where ilfc(^o) is the curvature contribution to Qtot today and Co ~ 3 is a cosmological- 
parameter dependent factor (see Appendix D). Observational bounds on curvature give 

pLS < 1. 

As mentioned above, perhaps the most interesting bubbles are those whose future light- 
cones intersect part of the reheating surface. The number of such bubbles can be obtained 
by integrating the distribution (4.18) over the range —pLS < < Pls- The total number of 
such bubbles is 

Nls = ^ [2PL5 + smh{27]Ls + 2pLs) - smh{2r]Ls - ^Pls)] (4.20) 
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Using Pls ^ i a-nd tjls ^ 1 this simplifies to 



Now using pls - Co\/^fe(^o) and tjls = log ^, we get 



Nls ~ ^PLse'"^'' (4.21) 



A^L5«^Co(^^) V^h^k) (4.22) 
Comparing this to the total number of bubbles, we have 

^ ~ KoV^kih) ~ l2y^fhM (4.23) 

valid for small ^/^l^JJo). 

For bubbles that intersect the last scattering surface, we can change variables from the 
distance £ to the angular size ip^s on the CMB sky. In other words, ipLS is the angular 
radius of the part of the last scattering surface affected by the coUision. Because the last 
scattering surface is small compared to the radius of curvature, we can approximate the 
spatial geometry inside the last scattering surface as flat. Also, the collision lightcone has 
grown large by last scattering, so it is nearly a straight line as it crosses the region inside 
the last scattering surface. Therefore we can use trigonometry to obtain 

i 

^ cos ijjLS (4.24) 

Pls 

Using the same approximations as above, the distribution in angles is 

d7V~^(^^^ y/n^)d{cosi;Ls)d'fl2 . (4.25) 

The distribution of angular sizes is rather featureless, 

dN oc d{cosipLs) ■ (4.26) 



5 Discussion 

Our primary goal in this paper was to re-consider the probability of observing a cosmic bubble 
collision in a model like the string theory landscape. We have found that incorporating the 
effects of a realistic cosmology inside the observation bubble has a dramatic effect on the 
result, due to the appearance of factors of the dimcnsionlcss ratio of the energy scale of the 
false vacuum to the energy scale of inflation inside the bubble. When this ratio is large 
(as it must be absent signiflcant flne-tuning) the number of bubble collisions in the past 
lightcone of an observer in a universe like ours is enhanced, ~ 7(if//ifi)^. Physically this 
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ratio arises for a simple reason: if the observation bubble inflates with Hubble constant Hi, 
late-time observers like us can see a region of the domain wall with surface area H~^, which 
is large in false vacuum Hubble units. As a result, we can see many Hubble volumes of the 
false vacuum in which collision bubbles may have nucleated. 

In addition to the total number we computed the collision distribution as a function 
of proper distance between nucleation points in the false vacuum, angular location on the 
observer's sky, proper time when the collision first enters the observers past lightcone, and 
angular size on the CMB sky (and on the wall of the observation bubble in Appendix B). The 
results are striking in several ways. The distribution, when projected onto the observation 
bubble wall at late times, is globally conformally invariant {i.e. invariant under the Lorentz 
transformations of the nucleation point) in the limit of small angular size. The anisotropy 
in the distribution, which is present for all but a special observer inside, turns out to be 
extremely small (because nearly all of the collisions visible to such an observer have very 
small size on the bubble wall, and the distribution in that limit is rotationally invariant). 
In particular, as we will see in Appendix C the anisotropic terms in the distribution are 
suppressed by a factor of (Hi/Hf)^. 

The distribution on the last scattering surface is of particular relevance for the potential 
observability of these events. Collisions with lightcones that bisect the part of the last 
scattering surface we can observe can create hot or cold disks on the CMB temperature 
map, and if these disks are of reasonable size and contrast they could be detectable [13]. We 
find that the number of such disks is N^s ~ ~ ^{Hf/HiYy/n^. Hence N^s scales 

as VHjfc = 1/d ~ Q~{N-N*)^ where N* is the number of cfolds of inflation needed to give the 
universe a radius curvature of order .1 in Hubble units (as required by current observational 
bounds) . 

As expected, long inflation inflates away the effects of bubble collisions by inflating away 
their number density — but since the experimental constraint on curvature is not very strong, 
the total number can still be quite large with no flne-tuning. However even if N^s > 1) one 
must also consider the effects of inflation on the magnitude of the signal from each collision. 
The analysis of [13] showed that a single collision with N — N* ^ 5 can have an easily 
observable effect without flne-tuning the parameters of the coUision bubble itself or our 
boost towards it (which is determined by X in our notation). If N^s is large or if some 
of the collisions are with bubble types that have particularly dramatic effects on the early 
universe, — A^* could be considerably larger than this and still produce observable effects 
on the CMB temperature map. 

The results of this analysis are interesting for a number of reasons. For one, an observation 
of the effects of a coUision of this type would be, if not a direct conflrmation of the string 
theory landscape, at least very dramatic evidence in its favor. The observational signatures of 
these collisions has been studied in detail in the context of a collision between two bubbles 
in a number of recent papers [12, 6, 13, 7], and if they are absent in the data this could 
constitute an interesting constraint on string models. 

Our flndings open a large array of directions for future research. One important issue 
is to make more explicit the connection between the thin-wall Coleman de-Luccia bubbles 
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considered here and the instanton transitions in the string theory landscape itself. Another 

is to further study the effects of the collisions on inflation and the CMB sky, and to in- 
vestigate additional observational opportunities like the effects on large scale structure or 
CMB polarization. A particularly interesting possibility along those lines is cosmic 21 cm 
data, which has tremendous potential as a probe of high-scale physics in the early universe 
[16, 17, 18]. It would be very useful to run a numerical simulation to verify the validity 
our assumptions about the domain wall dynamics and their effects of the reheating surface. 
There is an intriguing connection to conformal field theory, which will be explored in [19]. 
Finally, there remains the larger issue of the probability measure in the full string theory 
landscape. 

Appendices: In Appendix A, we address the issue of bubbles nucleating inside bubbles 
with a rate different than that of the parent false vacuum. We find that even in the extreme 
case that all bubbles are absolutely stable (so that no nucleations occur inside them), our 
distribution is unchanged up to factors of order jln^Hf /Hi), which is generically small. 

In Appendix B we compute the distribution of collisions on the asymptotic boundary 
of the observation bubble using a slightly different technique than that in the bulk of the 
paper. The results are consistent up to the expected accuracy, and we make the connection 
to the results in [4] more explicit. 

In Appendix C we quantify in detail the level of anisotropy in the collision distribution. 

In Appendix D we derive an expression for the conformal time including the effects of 
early-time radiation domination and late-time dark energy domination. 
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A Proper treatment of overlapping bubbles 

The analyses in the paper did not take into account the effects of the interactions of bubble 
nucleations with each other. After a collision bubble has nucleated, its future lightcone is no 
longer available for the same type of nucleation event. The vacuum inside the collision bubble 
may be completely stable; in any case it will certainly have a different set of decay rates than 
that of the parent false vacuum. As a simple test to see how significant this effect can be 
for our analysis we will assume the bubbles are absolutely stable, and therefore eliminate all 
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bubbles from our distribution that would have nucleated inside the future lightcone of other 
bubbles. 

Given two points A and B in de Sitter space, we need a formula for when A and B are 
spacelike separated. In terms of the embedding coordinates, the condition is 

V^.X'^X^s = XaXb + RaRb cos{5e) - TaTb < 1 (A.l) 

where 59 \s the angular separation between the centers of the two bubbles. Changing to our 
coordinates, the condition is 

cosh(rA — tb) < cosh(A'^ — Xb) + cosh ta cosh tb(1 — cos{59)) (A. 2) 

To take this restriction into accoTint, one could distribute bubbles according to (2.9) and 
then eliminate any bubble that is in the future lightcone of another bubble. However, we 
will show that the expected number of bubbles which must be eliminated is generically much 
less than one, so that this correction can be neglected. 

To show this, we consider the S2 defined by the intersection of the observer's backward 
lightcone with the domain wall of the observation bubble. In coordinates, this S2 is defined 
by ^7 = —00, p = 00, r] + p = r]Q. Since this S2 is on the lightcone of the observation bubble, 
it can also be defined in the coordinates that cover the exterior of the bubble as the surface 

X = —00, T = 00, X + T = Tjo. 

A given coUision bubble affects a solid angle Q on this 5*2. We can use the formula (A. 2) 
to find that the sohd angle inside the future lightcone of a bubble with coordinates {X, r]^) 
is 

n = 27r^^^- (A.3) 

cosh(?7^ — X) 

while the angular location of the bubble is randomly distributed. 

A sufficient (but not necessary) condition for two nucleation events to be spacelike sepa- 
rated is that the solid angle on the observation bubble's wall affected by one is not completely 
contained in the solid angle affected by the other. Therefore, if the expected total solid angle 
affected by all of the bubbles is much less than iir, then it is unlikely that any bubble is in 
the future lightcone of another. The expected total sohd angle is 

^47r^ V, ^'i^J_J^J^'' eosh^A" 2 cosh(r;„ - A") H 3 + ""^^ 

(A.4) 

As we have seen, for a realistic cosmology rio ~ In 7/. Since 7 is exponentially small, 7770 
will generally be much less than 1. 

B Collision distribution on the asymptotic sphere 

In this appendix we will compute the bubble collision distribution as a function of the angular 
radius on the asymptotic boundary sphere — which is the boundary of any open spatial slice 
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inside the observation bubble, or equivalently the observation bubble's lightcone after infinite 
time — taken up by collisions. This in general is not the same as the observed angular size 
of a bubble as seen by a finite time observer. In fact, not all the bubbles that chop off 
a portion of the boundary are going to be visible to any given observer even after infinite 
time unless Am = 0. Hence, the distribution we are going to compute here is not a directly 
observable quantity. Nevertheless we find it a useful concept, as it characterizes the statistics 
of collision events unambiguously, and independently of the observer and of the Hubble rate 
inside our bubble. Indeed, precisely because we are not keeping track of the propagation of 
the collision signal inside our bubble, the distribution we get cannot depend on the geometry 
of our bubble, but only on the geometry of its wall — which we approximate as the future 
light-cone of the nucleation event — and of the outside. As a result, this distribution will 
be applicable to any cosmology inside. (For the distribution actually seen by a realistic 
observer, see the analysis of Section 4.) 

For simplicity we set Hf = 1. The dependence on Hf can be recovered at any stage of 
our computation just by dimensional analysis. In fiat FRW coordinates the metric outside 
our bubble is 

It is actually more convenient to use the {a,F) coordinatization defined by [4], 

a = e* , F = r - e"* , (B.2) 

for which the metric reads 




Figure 4: The (a, F) coordinate system described in the text. Our bubble nucleates at the origin 
{t = 0,r = 0). 
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a is just the scale factor, while the constant-F hypersurfaces are past directed light-cones, 
with tip at r = (see fig. 4). This chart can in fact cover all of deSitter space, if we let a and 
F take any real values, with the constraint aF > — 1 (r > 0). However we will be interested 
in a smaller range of a and F. First, the surface where initial conditions are given (which we 
take following [4]) corresponds to a = 0. Then, our bubble's nucleation point, t = 0,r = 0, 
corresponds to a — 1,F = —1. Therefore, its past light-cone is the surface F — —1, with 
< a < 1. On the other hand, our bubble's wall — the future-directed light-cone with tip 
at (t = 0,r = 0) — is the surface a = 2/(1 — F), with — 1 < F < 1. The region outside 
both these light-cones, where other nucleation events can take place and eventually lead to 
collisions with our bubble is 

-1<F<1, 0<a< . (B.4) 

1 — r 

Notice in particular that F < 1 is necessary for the bubble nucleated at (a, F) to eventually 
collide with ours. 

Consider then a bubble nucleated at (a, F). Its angular size t/j on the boundary sphere 
is given by [4] 

1 g + 2F + aF^ 

Li = cosip= . (B.5 

^ ^ 2 1 + aF ^ ^ 

Since in the end we are interested in the distribution as a function of n, we can change 
variables and trade the nucleation a for /i, 

^ "-^ (B.6) 



1 + F2 - 2F/i ■ 

This way each collision event will be labeled by its angular size ^ and by F, which is 'con- 
served' between nucleation and collision — because the bubble wall follows a null trajectory 
with F = const. The allowed ranges for these variables are 

-l</i<l, -KF <n. (B.7) 

Indeed eq. (B.5) is an increasing function of a for fixed F. This means that the maximum 
and minimum of /i are attained at the extrema of a, for a given F. This yields the ranges 
above. The infinitesimal four-volume element corresponding to the infinitesimal intervals 
d/i, dF is 

dV^ = 47r(l + aFf da dF = Stt (1 - F ) 

^ ' (l + F2-2F/i)4 ^ ' ^ ' 

where we used eq. (B.6) for the change of variable a — > /x, and the iir comes from the solid 
angle integral. Then, the four- volume where a bubble of size between /i and fj. + dfj, can 
nucleate is just the integral of this along F, from —1 to /j, (see eq. (B.7)), 

^= /Vj^^4vr (2-;.) 

dn J_i diidF 3 (1 - /x)2 ■ ^ ' ' 
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A bubble of size /i takes out a solid angle 27r(l — //) on the boundary sphere. The sohd angle 
that is free of bubbles with angular sizes from /i — —1 to /i obeys the equation 

^ = -2.(1-/.). 7-^^, (B.IO) 

where 7 is the nucleation rate, and we included the factor Qf/An to account for the fact that 
at any given fj, only such a fraction of the solid angle is still available for new nucleations 
and collisions (in other words, we are counting regions in which disks overlap only once, 
no matter how many such overlaps there may be) . The solution with boundary condition 
= Att is 

^^/(/^)=47r^ e-^^(^+'') (B.ll) 



or, in differential form. 



dfj, 3 1 — // 



g-37r7(l+Mj (;b.12) 



As a check notice that if we consider all bubble sizes, from /i — —1 to /i — 1, all the solid 
angle but a set of measure zero is eaten up by collisions, 

Q/(l) = , (B.13) 

as expected. Another quantity of interest is the expected number of bubbles with angular 
size between /j, and 11 + d/i, 



dN{n) _ 7 rffi/(/i) 

dfi 27r(l — fi) dfi 

Integrating this we get total number of bubbles down to angular size /i 

47r7 / 3 ^ ^ 



(B.14) 



Nil,) = 



4e^7r7 

27r7 



7(^,-^(1-//),-^) 
7(- 1 + ^,-^(1-/.),-^) 



(B.15) 



where 7(0, zi, Z2) is the generalized incomplete gamma function, 7(0, zi, Z2) = J^^ t'^~^e''^dt. 
Perhaps more interesting is the asymptotic behavior of n{ii) at small angular sizes ^0 <^ 1. 
Expanding /i — cos ip we get 



which gives us the fractal dimension of our ensemble, d = 2 — in agreement with GGV's 
result [4]. 



22 



Although we did not keep track of the bubble positions on the boundary sphere, the 
distribution we have computed is obviously isotropic, because we have computed it in a 
frame where the initial-condition surface (a = 0) is invariant under rotations. This is the 
distribution that could in principle be seen by the central observer after infinite time, if she 
lives in a bubble of Minkowski vacuum and can see all the way to her bubble's domain wall. 
To compute the boundary distribution relevant for a different observer, we have to "boost" 
eq. (B.ll). 



B.l Trasformation properties under boosts 

We are thus led to consider how a deSitter boost acts on the boundary sphere. Recall that 
deSitter boosts about our bubble's nucleation point (henceforth "the origin") act as spatial 
translations on any given hyperbolic FRW slice inside the bubble. That is why different 
comoving FRW observers are mapped into each other by boosts about the origin. However 
since we are interested just in the boundary sphere, we can use the fact that the isometrics 
of three-dimensional hyperbolic space acts as the two-dimensional conformal group on the 
boundary. In particular, a translation (= deSitter boost) along z acts as a dilation about 
the south pole and a contraction about the north pole. More precisely, consider the metric 
on the sphere 

it is a conformally flat metric, as made evident by the change of variable 



1 + cos 6* „ 2/9 „, 

p=— — ^ sine=— ^, B.18 
sm fc' 1 + 



which is the standard stereographic projection, mapping the sphere to the plane. The North 
Pole (NP) ^ = gets projected to infinity, whereas the South Pole (SP) 6* = tt gets mapped 
to the origin p = 0. Upon the stereographic projection the metric becomes 

rf^^^=^Y^(rf/ + /V). (B.19) 

Then, dilations about the origin, p — > Ap, leave the metric explicitly conformally flat. These 
are the trasformations that correspond to translations along z in the three-dimensional hy- 
perbolic interior of the sphere. In terms of Q and <p they act as 

2A 

sin B sin ff = sin ^^^^^^^ . ^ ^ f' = f ■ (B^20) 

The action of this transformation is schematically depicted in fig. 5. 

Bubble collisions chop off circles on the boundary sphere. Interestingly enough, circles 
get mapped into circles by our conformal transformation eq. (B.20). To see this, recall 
that the stereographic projection (B.18) maps circles on the sphere to circles on the plane. 
Combining this with a dilation on the plane p — > Xp, which clearly maps circles into circles. 
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Figure 5: Left: A boost along z acts on the boundary sphere as a dilation about the south pole 
and a contraction about the north pole. Right: A circle on the boundary sphere is parameterized 
by its center's coordinates {6, (p) and its angular radius ijj. 



and with an inverse stereographic projection, we get that eq. (B.20) maps circles to circles. 
So, given a circle we just have to understand what happens to its center and its radius upon 
the transformation. Notice that the new center's 9 will not be simply the transformed 9 of 
the old center according to eq. (B.20). Rather, given a bubble with center at 9 and angular 
radius (see fig. 5), to find the new center's position we have to transform the two extrema 
of the bubble along theta, (^ ± ■?/'), and take their half-sum, 

ff . '^ + ^')';"'-^')' . (B.21) 
Likewise, the new angular radius is given by 

, (9 + i))' ~ (9 - 

i)' = ^—Jll ^ ZL (B.22) 

To simplify the algebra, we can consider two limiting cases: that of an infinitesimal 
transformation, A = 1 + e, and that of an infinite one, A ^ 1. In the former case we have an 
infinitesimal, 6'-dependent shift on 6', 

59 = 9' - 9 ^ -ei^m9 . (B.23) 

For a bubble centered at 9 with angular radius ip we get 

~ — e sin cos , ■?/'' ^ — e sin cos ^ . (B.24) 

If we restrict to very small bubbles, -^1, these transformation laws further simplify to 

^-^-esin^, ^' ~ ^/^(l -ecos^) . (B.25) 
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The number of bubbles per solid angle down to angular size ^ is given by eq. (B.16) 

dN —^d{cose) , (B.26) 

■0^ 3 

where the d{cos 9) comes from the isotropy of our distribution. The number dn is invari- 
ant under our boost, whereas 9 and transform as in eq. (B.25). Therefore the boosted 
distribution reads 

dN — ^ (l + e^f2cos^')d(cos^') ■ (B.27) 

7/y2 3 

The boost introduces a dipolar anisotropy in the bubble-collision distribution. 

More interesting for us is the opposite limit, that of infinite boost A — >^ oo. To figure out 
its action, it is best to work directly on the stereographic plane. Given a circle on the sphere 
with center at 9 and angular size ■0, this corresponds to a circle on the plane with center at 
Pc and radius p;,, 

Pe=i(p+ + P-), (B.28) 

where p± are the projections of ^± = 9 ± tjj, according to eq. (B.18). After some some 
straightforward algebra we get 



Pc = , Pb = , B.29) 

p — X p — X 

where we have defined 

X = cos 9 , ji — cos '4> . (B.30) 

The inverse mapping is 

^ ^ pI-pI-1 ^ pI-pI + i 

^{pl + pI + 1)' - ^pIpI ' ^{pl + pI + 1)' - ^pIpI ■ 

Now, the number of coUisions per unit x and p, is (see eqs. (B.14, B.12), and recall that 
before boosting the distribution is isotropic) 

2 

~ ^^f - ^(T^ (^^) dxdp , (B.32) 

which in terms of pc and reads 



dN ~ 7 



-1 \ 3777 



X ^ dpb dpc ■ (B.33) 

Pb 



The piece outside the brackets is manifestly invariant under dilations p \p, which are 
precisely the transformations we are interested in. We just have to figure out what happens 
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to the terms in brackets. Performing a very large boost corresponds to using a new radial 
coordinate 

p' = Ap , A > 1 . (B.34) 

Finite values of p' correspond to tiny values of p, which on the sphere correspond to very 
small displacements from the SP, x ~ —1 and p ~ 1. More precisely, from eq. (B.31) we 

have 

p ^ 1 - ^p;2 . (B.35) 
Plugging this into our distribution and dropping the primes we get 

dN^^ ^dpbdp, (B.36) 

A3'^^P63-3^7 

which projected back onto the sphere via eq. (B.29) reads 

dN ^ 4 ^cos^ _ 

As'^T (sinV')^"3^^ (cosV- - cos^)3'^^ 

To sec the relation to the distribution derived in the bulk of the paper (Eq. (2.9)), wc 
need to integrate the latter over X. If we do not impose any restriction on the observer from 
domain walls, the limits of integration are from (—00,00), giving (c.f. Eq. (2.9)) 

Stt 

dN = —7 cosh^ Tdrd cos 9 . (B.38) 

where we have used 1]^ — X = t. The relation between the coordinate r and the asymptotic 
angular radius ip is simply tanhr = cosV'; with a little algebra this becomes 

dN=^-^^d cos 9. (B.39) 

So far we have placed no restrictions on the direction of the boost. However if we require 
that the boost avoid bubble domain walls, the distribution Eq. (B.37) will be modified by 
a factor of 6(6' — ^), which enforces that no disks cover the point ^ = the boost was 
taken towards. This condition coincides precisely to the one imposed by the initial condition 
surface at infinite boost discussed in the text below Eq. (3.8), which is cos'^ + tanhr > 0, 
or in these coordinates 9 > ip. 

Expanding Eq. (B.37) in powers of 7 gives 



dip 
siv? ip 



dN r^^—^d cos 9 {1 + {-f]nS)), (B.40) 



in agreement with Eq. (B.39) up to factors of order 7 In 5. Here 5 is a UV cutoff applied 
both to the minimum size ipram ~ 5 of the disks and to the boost parameter Amax ~ 1/5- 

The cutoff on boosts is necessary because Eq. (B.37) is an average derived under the 
assumption that no bubbles nucleate inside other bubbles. Boosting the resulting distribution 



26 



infinitely leads to a singular result, since the average was taken before the boost and all but 
a set of measure zero of such boosts are in directions already "eaten" by a bubble domain 
wall. 

A cutoff on ip corresponds to a cutoff on the observer's conformal time rjo through X+t = 
r)o. After integrating over X this translates into — In'^min ~ ^o; or 'V'min Hi/Hj. As a quick 
check, according to Eq. (B.37) the fraction of the asymptotic sky covered by disks is then 
7 J ip'^dip/ sin'^'f/' ~ — 7ln'?/'inin ~ 7770, in agreement with Eq. (A. 4). Hence, the expansion 
in 7 shows that the distributions (B.37) and (B.39) agree as expected up to corrections of 
order ■yrjo ~ jlia^Hf /Hi). 



C Legendre moments of the collision distribution 

In this appendix we derive the Legendre moments of the bubble collision direction distribu- 
tion, and compare this with the distribution found in Ref. [4]. We focus on the effects of 
inflation with Hi ^ Hf inside the bubble, but these results can be straightforwardly gener- 
alized to include arbitrary cosmological evolution inside the bubble. Using the restriction of 
Eq. (3.6) we can show that for an internal de Sitter space with Hi ^ Hf the four-volume 
available for the nucleation of bubbles in the ^ — > 00, infinite boost, limit is 



dtdn sni [Hi cosh^ (¥) - sinh^ (^)] 

nlsux^Hit) 



X 



n1{l - cose) cosh^ (^) + (1 + cos^) sinh^ (^) 



— tanh 



Hit 



(C.l) 



where we have adopted the notation Tij = Hi/ Hf for the dimensionless ratio of Hubble rates. 
In the limit Tii ^ 1 this reduces to Eq. (45) of Ref. [4] . 
We define 

jjL — cos 9 , 



and write 



where 



dVA 



a 



a — 



dtdn {p - At)3 
Hf sinh^ (Hit) 



-A 



P 



and 



A 



3n|[njcosh^^)-smh^^f)] 

Hfcosh^ (^)+sinh^(^) 
?ifcosh2(^)-sinh2(2|*) 

tanh^ (^) 
37if [7ifcosh2(^)-sinh^ (^)]' 



(C.2) 
(C.3) 

(C.4) 
(C.5) 

(C.6) 
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as 



We expand Eq. C.3 in orthonormal Legendre polynomials that satisfy 

J d/lPMPmilJ') = Smn 
dV4 



dtdQ 



n=0 



The generating function for the Pn(/^) is 



n=0 



2n + 1 



so that 



(C.7) 
(C.8) 

(C.9) 
(C.IO) 



Prom Z{^, t) we can define a generating function for the Legendre moments Vn as 

^ I I \/ A I ^1 \ I I. 'II. I 

(C.ll) 
-2A, (C.12) 



a 



zf-^'p-T + 34«'arcoth[Z;,(/3 - u)] 



where = Z{(3, u), such that 



u=0 



The first few moments are 



^0- (/32_i)2-^ 



^^1 



V2 



V5a{5p - 3/33 + 3(/52 - l)2arcoth/3) 



V2(/32 - 1)2 



in the limit Hd ^ 1 we have 



^^0 



2\/2 



" 2 1 



(C.13) 



(C.14) 



(C.15) 



(C.16) 



(C.17) 
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Vl 



2^6 



" 1 



-Hit 



+ 0(e-2W>*) 



(C.18) 



V2 



2V10 



-1 + 87Yf + Unf In 7^1 - 87if + 



(C.19) 



from which we can see the dipole and higher multipoles are suppressed relative to the 
monopole by a factor of in the hmit that TYj is a small parameter, and the distribu- 
tion is nearly isotropic. 

To leading order in the Tii expansion we can write for an arbitrary multiple moment 



vn = (-1)"+^ ^%^"^^ [1 - nl^n + 2)(n - 1) + e'^^* + O(e-2«'0 (C.20) 

from which we see that all multipole moments are present, albeit suppressed by a factor of 
Tif, at late times when Hi ^ 1 rather than the late-time monopole- plus-dipole distribution 
found in Ref. [4] for the Hi = 1 case. 

The late time distribution takes the exact form 



cIVa 



1 - 



8H^ 



dtdfl 3HI{1 - Hi) 
which in the limit 7i j — > 1 is simply 

dtdVt ' 



[l + cos^ + 7if(l-cos^)]' 



(C.21) 



2(l + cose)e-* + 0(e-2*) 



(C.22) 



in agreement with Ref. [4]. 

To leading order in Hi the late-time distribution is 



dV. 



dtdft mf 



1 + H^ + Hj+ \1 



[1 + cos^ + 7^1(1- cos ^)]' 



+ 0{H^) 



(C.23) 



and so the anisotropic part of the distribution is suppressed by Hf except in the region 
n — 6 < €, where e Hi/Hj. Only a fraction of order / Hj of the collisions are centered 
within this region. 



D Computation of the conformal time with a 7^ ^fc 7^ 7^ 

In this appendix we derive an expression for the conformal time taking into account early time 
radiation domination and late time dark-energy domination (modeled here as a cosmological 
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constant) of our Universe. In this paper we use a convention where the scale factor has units 

of length and the curvature A; = —1 is a dimensionless quantity. We will not set A; = —1 in 
this derivation to allow easy comparison to the alternative convention often used where the 
scale factor is dimensionless and k has units of inverse length squared. For a scale factor 
a{t), where t is cosmological time, the dimensionless conformal time today since some early 
time ti is 

/•*(' dt' 

Ari = rio-rie = V-k / —jr , (D.l) 

which can be expressed in terms of the Hubble rate H as 

P da 

Arj^V^ -—. (D.2) 



. am 

We introduce the variables a — a/a^ where a^ is the scale factor today, and H — H/Hq 
where Hq is the Hubble rate today. In terms of these variables the Hubble rate is 



H{a) = ^Vtr a-4 + VL^ d-^ + Vtk a-^ + (D.3) 

and so 

/(a) = ^i — ^ ^ -. (D.4) 

a?H{d) v^r + S + fife + a* 

The conformal time is then 

^^-^ f ^«/(«) -V^k f daf{a) . (D.5) 

Oo-no Jai Jai 

where we have used the definition of the curvature density Qfe = —k/{a'^Hl). We write 

C,Q ^ CQ{Vtr,Vtm,^k,^h) ^ I daf{a) (D.6) 

J hi 

Typically ai is the scale factor at reheating and we can safely take czj ~ 0. We set ^ in 
what follows. To make further progress we split the integral about Og = {VLr/^m)-, the scale 
factor at matter-radiation equality as 

Co ~ r d~af{a) + / d~af{a (D.7) 

JQ Jae 

As it starts deep in the radiation dominated epoch and ends at matter-radiation equality 
the first part of the integral can be very accurately approximated in the limit fl/i and 
f2fe — > as 
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ddf{a) ~ ( V2 - 1) 



The second part of the integral can be solved as a series expansion about 
write 



(D.8) 
ilr = 0. We 



1 gim+n) j 



mini dVL'^dVt't 

m=0 n=0 ' 



Qfe=0 

Qr=0 



E 



1 d''f 



n=0 ' 



df 



n,.=o 



+ 



(D.9) 



where in second line we have only included the terms that contribute at leading order. We 
perform the integral term by term to find 



„1 oo 

d~af{d) = 

J 5,^ , n 



where 



n=0 



2ra+l 4ra+l 4n+4 



(-l)"+i(2n- 1)!! 1 
2"(4n + l)n! 



and 



3fi 



3/2 



Ha 



1 



1 



;4n+l 



2^1 



2n+l 4n+l 4n+4 



-1/2 



- 1 + 



-1/2 



(D.IO) 



(D.ll) 



(D.12) 



Combining these expressions, expanding to leading order in the cosmological parameters, 
and substituting cie = flr/^m we find that 



nm 



1 



1/3^1/6 



1 - 



80a 



(D.13) 



to a very good approximation (higher order terms contribute < 0.1% near current cos- 
mological parameters). Note that while Eq. (D.IO) is valid for any value of Qa including 
f^A — ^ 0, Eq. (D.13) is a good approximation only in a Universe, like our own, in which the 
cosmological constant is already dominating the expansion. 

Current cosmological observations [20] admit 0^ ~ 0.01. For instance, recent cosmo- 
logical constraints are compatible with Qa = 0.72 ± 0.03 and Jl^ = 4.12 x 10~^//i^ ~ 
(8.3 ± 0.6) X 10-5 near Vtk ~ 0.005 or Qa ^ 0.71 ± 0.01 and Vtr ~ (7.9 ± 0.2) x 10"^ 
near ~ 0.01. These translate to a typical ranges Co — 3.35lg }4 near VLj^ — 0.005 and 
Co — 3.33^*^0 05 near fij. ~ 0.01 — more than 60% larger than the matter-dominated estimate 
(C — 2) discussed in Section 4.1 of this paper. 
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